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Abstract 

A special series is introduced in this paper to yield solution of the first-order linear 
vector differential equation. It is proved that if the differential equation satisfied by 
the first term of this series can be solved exactly, then other terms can be determined 
by the method of variation of parameters. We point out that the special series will be 
the solution of the first-order linear vector differential equation if the infinite special 
series converges. An illustrative example has been given to outline the procedure of 
our method. 
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1 Introduction 


The content of linear systems constitutes a large and very important part 
in the theory of ordinary differential equations. From the early days of ordi¬ 
nary differential equations the subject of linear systems has been an area of 
great theoretical research and practical applications, and it continues to be so 
today. In physics, many dynamical systems or material helds can be treated 
as linear systems. For instance, the evolution equation of a spin-| fermion in 
Robertson-Walker space-time deduced from the covariant Dirac equation is a 
linear system [1]. In this paper we shall restrict our attention to linear systems 
of n differential equations in n unknown functions only. 
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Almost in all books on ordinary differential equations, the authors always 
present a symbolic operator method for solving linear systems with constant 
coefficients. After introducing the so-called differential operators, the standard 
procedure for solving linear systems with constant coefficients will be described 
in any such book [2]. Therefore, we shall only focus on the linear systems whose 
coefficients are functions. 

When the system of linear differential equations is of the type that have 
variable coefficients, the method of solving it depends upon the concrete form 
of the coefficients and the type of differential equations. One can try to solve 
this kind of linear systems by acquiring a higher order differential equation 
which can be solved [3,4], or to solve them approximately by various meth¬ 
ods such as numerical methods, series methods, and graphical methods. In 
this paper We shall introduce a relatively general and systematic method for 
solving linear systems with variable coefficients. 

The present paper is devoted to a method for obtaining solution of the 
first-order linear vector differential equation in special series form. We shall 
decompose the linear vector differential equation into a group of correlated 
vector differential equations. After this, we shall proceed to explain the method 
of special series. We shall then discuss the convergence of infinite special series, 
and point out the solution. Finally, we shall give an illustrative example to 
clarify the procedure of special-series method. 


2 Decomposing the linear vector differential equation 


We consider the normal form of linear system of n first-order differential 
equations in n unknown functions j/i, |/ 2 ,.. .,?/n- It is well known that this 
system is of the form 

y = A(x)y+ F(x) , (1) 


where we have introduced the notation y' = ^, the vectors y and F(x) are 
defined respectively by 





^ h{x) ^ 


?/2 


f2{x) 

y = 


and F(a;) = 



\yn) 




( 2 ) 
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and the matrix function A(a:) being defined by 


A(a;) = 


^aii(x) ai2{x) ■ • ■ ain{x) ^ 

021 (x) 022(2:) ■ • ■ a2n{x) 


y 0^1(2:) 0^2(2:) ■ ■ ■ CLnnix^ j 


( 3 ) 


We shall assume that all of the functions defined by fi{x), i = 1,2, 
and aij{x), i = 1, 2,..., n, j = 1, 2,..., n, are continuous on a real interval 
[a, b]. If all fi{x) = 0, z = 1, 2,..., n, for all x, then the system (1) is called 
homogeneous. Otherwise, the system is called nonhomogeneous. 

Hence the normal form of linear system (1) is a nonhomogeneous linear 
vector differential equation, and the corresponding homogeneous linear vector 
differential equation reads 

y = A(x)y . (4) 


Mathematician has developed the method of variation of parameters for find¬ 
ing a solution of the equation (1), assuming that one has known a fundamental 
matrix of the corresponding equation (4). Therefore, in this paper, we will fo¬ 
cus on explaining the method of special series for finding the fundamental 
matrix of the homogeneous linear vector differential equation (4). 

We shall assume that A(x) is a nonsingular n x n matrix function. Then 
one can divide A(x) into two nxn nonsingular matrices H(x) and Z(x). It is 
stressed that the elements of matrix function Z(x) must be hnite on the real 
interval [a, h]. That is, 

A(x) = H(x) + Z(x) . (5) 


Here H(x) can also be a constant matrix. For brevity, we use the common 
sign H(x) to denote both a matrix function defined on all real x G [a, b] and 
a constant matrix. Therefore, the equation (4) becomes 

y'= H(x)y-F Z(x)y . (6) 
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To find the fundamental matrix of the above equation, we shall assume that 
the vector function y is a special series written as follows 

CXD 

y = ^yi- (7) 

i=0 

And we shall assume that the first term yo of this series satisfies 

yo = H(x)yo . (8) 

Under the assumption that y is the sum of the first two terms, inserting 
y = yo + yi into the equation (6) yields 

yo + y'l = H(x)(yo + yi) + Z(x)(yo + yi) . (9) 

With the help of the equation (8), and ignoring the term Z(a:)yi in the equa¬ 
tion (9), we obtain 

yi = H(x)yi + Z(x)yo . (10) 

Under the assumption that y = yo + yi + y 2 , inserting yo + yi + y 2 into the 
equation (6) yields 

yi) + yi + yi = H(a:)(yo + yi + y 2 ) + Z(a:)(yo + yi + y 2 ) • (H) 

With the help of the equations (8), (10), and ignoring the term Z(x)y 2 in the 
above equation, we obtain 

yi = H(x)y 2 + Z(x)yi . (12) 

Under the assumption that y = jyjLoyj, inserting jyjLoyj into the equation 
(6) yields 

m m m 

y'j = h(2:) yj + yj ■ (i3) 

j=0 j=0 j=0 

Using the same method as before, and ignoring the term Z(x)ym in the above 
equation, we acquire 

y'^ = H(a;)ym Z(x)ym-i • (14) 
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Through the foregoing logic chain, we can draw the following conclusion: in 
the case of homogeneous linear vector differential equation (4), after expanding 
the vector of functions y into an inhnite series, at the same time, dividing the 
coefficient matrix into two nonsingular matrices, we can obtain a recursion 
series, each term of which satisfies 


I y' = B.{x)yj j = 0 

\ y'j = H(a:)yj + Z{x)yj-i j = 1, 2, 3,... 


In the case of nonhomogeneous linear vector differential equation (1), we can 
similarly dehne a special series, each term of which satishes 


I y' = B.{x)yj + F{x) j = 0 

\ y'j = H(a:)yj + Z(x)yj_i j = 1, 2, 3,... 


(16) 


From the theoretical point of view, we can say, the above equations demon¬ 
strate that one can obtain every terms in the recursion series y = y* 
and only if the first term can be solved explicitly. 


3 The solution of each term of special series 


The theory of linear systems has proven that if the vector functions <Fi, $ 2 , 
... ,$j are i solutions of (8) and Ci, C 2 , ■ ■ ■■, Ci are i numbers, then the vector 
function 


^ , 
k=l 


(17) 


is also a solution of (8). 

We shall be concerned with n vector functions, and we shall use the following 
common notation for the n vector functions in the following discussion. We 
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let $ 1 , $ 2 , • • • ,‘^’n be the n vector functions defined respectively by 



f^n] 


A.] 


\ 

*hi(x) = 

021 

, *^>2(2;) = 

022 

, . .., <l>„(x) = 

02n 


^ 0nl y 


y4>n2 y 


^ 0nn J 


( 18 ) 


The n X n determinant 


011 012 • • • 01n 
021 022 • • • 02n 

0nl 0n2 • • • 0nn 


(19) 


is called the Wronskian of the n vector functions $ 1 , $ 2 , • • • dehned by 
(18). We will denote it by hh($i, $ 2 , • • •, d)„) and its value at the point x by 
W($i,$2,...,$n)(x). 


Let the vector functions $ 1 , $ 2 , • • • ,^n dehned by (18) be n solutions of the 
homogeneous linear vector differential equation (8) on the real interval [a,b]. 
These n solutions $ 1 , $ 2 , ... ,$n of (8) are linearly independent on [a,b] if 
and only if 

ty($i,$2,...,<h„)(x)7^o (20) 


for all X G [a,6]. A set of n linearly independent solutions of (8) is called a 
fundamental set of solutions of (8). If the vector functions $ 1 , $ 2 , ... ,*hn 
dehned by (18) make up a fundamental set of solutions of (8), then the n x n 
square matrix 


M(x) = 


'^01l(x) 0l2(x) . . . 0ln(x) ^ 

02i(a;) 022(a:) • • • 02n(a:) 


( 21 ) 


0n2(^) • • • 0 nn(^) J 
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is a fundamental matrix of (8). Furthermore, let $ 0 ( 2 ;) be an arbitrary solution 
of (8) on the real interval [a, b]. Then there exists a suitable constant vector* 

c = (ci,C 2 ,... ,Cn)^ (22) 


such that 

$o(x) = M(x)c 


(23) 


on [a,b]. 

Inserting the solution (23) of yo into the equation (10) yields 

y; = H(x)yi + Z(x)M(x)c . (24) 

We shall now proceed to obtain a solution of above equation by variation of 
parameters. It has been proven [2] that <l)i dehned by 

X 

<li(x) = M{x) J M-\t)Z{t)M{t)cdt , (25) 

XQ 


where Xq G [a, 5], is a solution of the nonhomogeneous linear vector differential 
equation (24) on [a, b]. 

We have shown that $1 is a solution of differential equation (10). Directly 
inserting this solution into the differential equation (12) yields 

y 2 = H(x)y 2 + Z(x)$i(x) . (26) 


We shall also proceed to obtain a solution of above equation by variation of 
parameters. It has been proven that $2 dehned by 

X t2 

l> 2 (x) = M(x) J M-\t 2 )Z{t 2 )M{t 2 ) J M-\ti)Z{U)M{ti)cdtidt 2 , (27) 

XQ XQ 


where Xq G [a, 6], is a solution of the nonhomogeneous linear vector differential 
equation (26) on [a,b]. 

* We denote the transpose of A by A^, where A being any matrix. 
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According to above argument, we can acquire $ 3 (x), $ 4 ( 2 :),... etc term by 
term. Now we shall assume that the solutions <ho(a^)i ^’ 1 ( 2 ;),..., $m-i( 2 ;) have 
been acquired. Then, from the recursion formula (14), we can simplify the 
differential equation satished by 

y'^ = H( 2 ;)ym + Z( 2 :)l>m-i( 2 ;) . (28) 


It is easy to prove that dehned by 


X tm 

= M. I J M-_‘ , 

XQ XQ 

t2 

• ■ • Mt2 J ■■■dtjn , 

XQ 


(29) 


where = M(t) and Xq G [a, b], is a solution of the nonhomogeneous linear 
vector differential equation (14) on [a,b]. 

Hence the above reasoning demonstrate that dehned by 


= mN M-‘^ ... 

XQ XQ 

^2 

J M~^Zt,Mt^cdtidt2 ■ ■ ■ dtj , j = 1, 2, 3,... 

XQ 


(30) 


where xq G [a, 6], is a solution of the nonhomogeneous linear vector differential 
equation (15) on [a,b]. 

It is well known that 4 * 0 ( 2 ^) dehned by 


$o(x) = M(x) J M-\t)F{t)dt , (31) 

XQ 


where Xq G [a, 6], is a solution of the nonhomogeneous linear vector diherential 
equation 


y'o = H(i)yo + F(i) 


(32) 


on the real interval [a,b]. A theorem [2] has been proven that an arbitrary 
solution $ 0 ( 3 ^) of fhe nonhomogeneous differential equation (32) is of the form 


$ 0 ( 2 ;) = *^* 0 ( 2 ;) + M(x)c 


(33) 


for a suitable choice of ci, C 2 ,..., c„. According to the same logic chain as 
what has been explained in the study of the homogeneous linear vector differ¬ 
ential equation (4), A solution of the nonhomogeneous linear vector differential 
equation (16) on [a, b] can easily be acquired, that is 



(34) 



where Xq G [a,b]. 


4 The solution of the first-order linear vector differential eqnation 

THEOREM the special series is the solution of the hrst-order linear 
vector differential equation if the inhnite special series converges on the real 
interval [a,b]. 

proof. Consider the partial sums of vectors [5] 

So = yo, 

(35) 

S 2 = yo + yi + y2, 


Then Si, the hh partial sum of vectors, is given by 


Si = yo + yi + yz H-H yi = H yt ■ 


(36) 


k=0 
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In the case of homogeneous linear vector differential equation (4), from the 
equations (15), we directly acquire 

S; = H(x)Sz + Z{x)Si_i , 1>1 . (37) 

In the case of nonhomogeneous linear vector differential equation (1), summing 
all the equations in (16) for the set 0 < j < /, we obtain 

S; = H(x)Sz + Z(x)Si_i + F(x) , 1>1. (38) 

Since the integer I is unrelated with variable x, the limit with respect to I and 
the differential with respect to x are commutative 

lim S( = (lim Si)' , (39) 

l—^OO >CO 

With the help of equation (39), seeking for the limit of the partial sum of 
vectors, for the homogeneous linear systems, from the equation (37), we obtain 

(lim S;)' = H(a:) lim Si + Z(x) lim S;_i , (40) 

>CO /—>CXD 1^00 

and for the nonhomogeneous linear systems, seeking for the limit of the partial 
sum of vectors in the equation (38) yields 

(lim S;)' = H(a:) lim Si + Z{x) lim S^-i + F(x) . (41) 

1^00 1^00 1^00 

Assuming that the inhnite special series converges on the real interval [a, b ], 
then 

lim Si = lim S;_i = S , (42) 

/—^■OO 1^00 

where S is a finite function on [a, b]. Therefore, inserting the above expression 
into the equation (40) yields 

S' = H(a:)S + Z(x)S = A(a:)S , (43) 


that is, the limit S defined by (15) and (42) satishes the homogeneous linear 
vector differential equation (4). Furthermore, inserting the formula (42) into 
the equation (41) yields 

S' = H(x)S + Z(x)S + F(x) = A(x)S + F{x) . (44) 
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The above equation demonstrates that the limit S defined by (16) and (42) 
satisfies the nonhomogeneous linear vector differential equation (1). Since the 
infinite series converges and has sum S if the sequence of partial sums con¬ 
verges to S, the special series converges to the solution of the first-order linear 
vector differential equation. 

Thus we have shown the validity of theorem. 

CONCLUSION If the special series converges, then: 

(1) an arbitrary solution 4>(a;) of the homogeneous linear vector differential 
equation (4) on [a, b] can be expressed as 

OO 

$(x) = ^l>,(x) , (45) 

j=0 


where ^j{x) is defined by (30); and 

(2) an arbitrary solution <I)(x) of the nonhomogeneous linear vector differ¬ 
ential equation (1) on [a,b] is of the form 

OO 

4(x) = i:4,(x) , (46) 


where 4’j(x) is defined by (34). 

Example 

Solve the system 

{y[ = (e-2^ - 36^^ + 2)yi + (e-^* - + 3)y2 

[ 1/2 = — l)yi + (3e^* — — 2)^2 


(47) 


on the real interval [0,1]. 

Clearly this is of the form (6), where 

H(x) = 


^2 3 ^ 

, Z(x) = 

^ - 3e2^ 

g-2x _ 9g2x \ 

V-1 -2y 



3g2x _ g-2x J 


( 48 ) 
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The corresponding homogeneous differential equation satished by yo is 


dyo 

dx 



( 49 ) 


we hnd that 


' 3e^ \ 

and 4 > 2 {x) = 

^ e " \ 





(50) 


constitute a fundamental set (pair of linearly independent solutions) of (49). 
Thus a fundamental matrix of (49) is given by 


M(a:) 



(61) 


From the fundamental matrix M(a:), we find that 


M-^(x) 



-3e^ 


(62) 


Thus the term M ^{x)Z{x)M{x) in the formula (30) becomes 


M-^(a;)Z(x)M(x) 


^0 2 \ 
v2 oj 


(53) 


We have indicated that an arbitrary solution of (49) can be presented as 


<ho(x) = M(x) 



(54) 
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where ci and C 2 are constants. Inserting the equation (53) into the formula 
(25) yields t 


l>i(a:) = M(a: 



( 66 ) 


One can acquire $ 2 ( 2 ^), ‘h 3 (x), ... from the formula (30) term by term, that is 

/ {2xf 


$ 2 ( 3 :) = M(a:) 


V^C2^ 


$ 3 ( 3 :) = M(3:) 


3! 


-C2 


(2a:)3 

VHrCi 


( 66 ) 


Obviously, the special series given by (54), (55), and (56) is convergent. By 
theorem, we know that 


/ 

$(x) = M(x) 

V 

= ^M(3:) 


I o I (2a;)2 {2xf , (2x)^ , \ 

Cl + 2xC2 + 2! ^ 3! C2 + 41 Cl + • ■ • | 

I o I (2a:)^ I ( 23 :)® I (2x)'* , I 

C 2 + 2xci + ^ 2! O + 3! Cl + ^ 4/ C 2 + • ■ • y 

^e^^(ci + C 2 ) + - C 2 ) \ 

ye^^(c2 + Cl) + e-^^(c2 - Cl) J 


is a solution of the linear system (47) for every real number x G [0,1]. Fur¬ 
thermore, we simplify the solution of (47) as follows 


$( 3 :) 


/ 

V 


(-e-3^ + 3e-^ + + 3e3^)ci + (e-^^ 

(g-3a: _ g-x _ gX _ _j_ ^_g-3a: 


— 3e -t- -I- 3 e^^)c2 

_1- g-x _ gX _ g3x^g^ 


where 2ci = ci and 2c2 = C 2 . Inserting the above vector into the linear vector 
differential equation (47) really produces two identities. 


An important fundamental property of a normal linear system (1) is its 
relationship to a single nth-order linear differential equation in one unknown 
function. The so-called normalized nth-order linear differential equation is of 
the form 

1 In this example, we have selected xq = 0 for simplicity. 
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It is well known that the single nth-order equation (57) can be transformed into 
a special case of the normal linear system (1) of n equations in n unknown 
functions. Hence the method of special series can be used for solving the 
normalized nth-order linear differential equation. 


5 Conclusions 


A special series has been introduced to yield the solution of the linear vector 
differential equation. A group of recursion linear vector differential equations 
satished by the terms of special series has been obtained and solved. We proved 
that the special series is the solution of the hrst-order linear vector differential 
equation if the inhnite special series converges. We have given an example to 
outline the procedure of our method. In principle, the method of special series 
can be used for solving any hrst-order linear vector differential equation if and 
only if the considered special series converges. 
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